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Several  specific  instances  of  the  diffraction  of  acoustic 
and  electromagnetic  waves  on  wavy  surfaces  of  various  types  are 
examined  on  the  basis  of  general  theory.  Here  the  finite  nature 
.of  the  dimensions  of  the  area  on  which  the  diffraction  occurs  is 
taken  into  account.  In  the  case  where  the  uneven  condition  is 
sinusoidal,  the  results  obtained  by  two  different  methods  are 
compared. 

In  work  [l](we  shall  designate  it  1  in  the  future)  an  approximate 
theory-  of  the  diffraction  of  acoustic  and  electromagnetic  waves  on  an 
uneven  surface,  which  had  uneven  areas  which  were  large  in  comparison 
to  the  length  of  the  wave,  on  the  assumption  that  the  surface  was  Infin¬ 
itely  extended,  but  that  the  uneven  conditions  were  periodic,  was  given. 
Ke  shall  consider  the  finite  nature  of  the  surface,  as  well  as  examine 
several  specific  examples,  below. 


1.  Oiffraction  on  an  Uneven  Surface  of  Limited  Dimensions 

Let  a  limited  sector  on  an  uneven  surface  participate  in  the 
scattering  of  waves.  This  may  take  place  both  as  a  consequence  of 
the  limited  nature  of  the  area  and  as  a  consequence  of  the  limited 
nature  of  its  illuminated  sector.  For  the  sake  of  simplicity,  we 
shall  examine  here  a  two-dimensional  problem,  i.e.,  we  shall  assume  that 
the  "active''  sector  of  the  uneven  surface  has  the  form  of  an  infinite 
ptrip  positioned  in  plane  xy  perpendicularly  to  the  wave  incidence  plane 
xz.  As  in  [1]  we  shall  describe  the  unevenness  as  a  function  depending, 
in  the  case  under  examination,  only  on  one  variable  A'.  Iv'e  shall  assume 
the  width  of  the  area  to  be  large  in  comparison  to  the  period  of  uneven¬ 
ness  Ax. 

For  an  example,  we  shall  examine  the  potential  of  the  acoustic 
field  of  a, scattered  wave,  which  in  the  case  of  an  infinite  area  is 
given  by  the  integral  (12, l)1  </'e  shr.l  introduce  as  an  integral  the 

1  In  the  case  of  an  electro:,..  _  n  -:ci <•  w&ve,  the  integral  (32,1)  is  equivalent 
to  it. 
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function  F(X) ,  which  characterizes  the  change  in  the  "illumination'* 
of  the  area  in  direction  X.  As  a  result,  if  one  considers  that  we 
are  examining  a  plane  problem  and  the  double  Fourier  series  is  reduced 
to  a  single  series,  we  shall  obtain 


+09 


F{xi  '^p  V  (**«  +  hy + Ml x 


'  x  exp  —  k.r,  +  mp)  X  +  i(uq  —  ky)  K)  dX  d  Ydkjlky.  ^  ^ 


If  the  area  is  illuminated  or  sonicated  with  a  sharply  limited 
beam  with  a  constant  intensity  throughout  its  cross-section,  or 
the  illumination  is  accomplished  by  an  infinite  plane  wave,  but  the 
area  itself  has  the  form  of  a  strip  beyond  the  limits  of  which  there 
is  no  reflection,  F(X)  will  be  equal  to  1  within  the_ strip  and  zero 
outside  of  it.  Generally,  the  illumination  function  F(X)  is  more  complex. 
Here  we  shall  consider  it  to  be  little  changing  over  the  wave  length  of 
unevenness  As. 

Integration  in'  Cl)  for  X  and  is  accomplished  in  the  same  manner 
as  in  §  2  and  3  (see  1)  for  an  infinite  area,  which  is  natural,  insofar 
as  the  area  is  assumed  to  be  infinite  indirection  y,  as  was  previously 
the  case.  As  a  result  we  shall  obtain 


<?”tK  2  \  B  (**•  h)  exp  {i  [kxx  +  ktz))  dkx  x 

x  j  F(X)exp  ■  ft*  +  mp)  X)  dX. 

—CO  * 


We  shall  write  the  expansion  of  the  function  F(X)  in  the  Fourier  integral: 


+® 

F(X)  =  \  (*)</«, 

—00 


(3) 


where  is  the  function  characterizing  the  spectral  expansion 
density  defined  by  the  equality: 

*(*)*•  &\e**FVQdX.  (4) 


2 


- - ■*-£*-***,•**  f^wC  ,-  '  " 


In  Fourier  integral  theory  a  relationship  is  known  which  is  analogous 
to  the  relatiorhip  o£  indeterminacy,  according  to  which  intervals  D  and 
Ax  of  the  changes  in  variables  X  and  x,  where  the  functions  F(X)  and  4>tx> 
have  values  notably  different  than  zero,  are  related  by  the  relationship: 


DtM  'v  2v 


(5) 


.  Since  we  are  suggesting  that  D  »  (many  periods  of  sinuosity  are 
being  illuminated),  the  latter  relationship  yields 


Ax  «  2ir/A  =  p  ,  (6) 

which  we  shall  utilize  in  the  future. 

Besides  that,  if  F(X)  is  a  function  having  a  single  maximum  in 
the  middle  of  the  area,  then  it  follows  from  (4)  that  4>(x)  will  be 
maximum  when  tc  =  0.  Comparing  (2)  and  (4),  we  find 

<?’=  2  —  +  mp)exp[l(k*x  +  kiz))dkjl. 

/«»<*— co  — <x>  ^  J 


From  here  it  is  evident  that  in  the  case  of  a  limited  area  there  are 
scattered  waves  with  arbitrary  kx,  i.e.,  waves  distributing  themselves 
in  all  directions.  However,  for  each  m  the  amplitudes  of  these  waves 
differ  notably  from  zero  only  within  the  limits  of 

kx  --  k°x  ■*  mp  +  Ax  , 

since  beyond  these  limits  the  function  -  kx  +  nip)  turns  out  to 

be  practically  equal  to  zero.  Since  according  to  (6)  Ax  «p,  it  consequently 
follows  that  the  scattered  waves  will  group  themselves  in  narrow  angular 
intervals  close  to  the  directions  set  by  the  equation 


kx  =  k°x  +  mpt  m  =?  0 ,  +_  1 ,  +_  2 ,  ...  . 

We  shall  examine  the  field  at  sufficiently  great  distances  from 
the  area,  par  .icularly  in  the  Fraunhofer  zone.  We  shall  make  the 
beginning  of  the  coordinates  and  the  center  of  the  area  coincident 
and  we  shall  define  the  angle  created  by  the  direction  to  the  point 
of  observation  and  the  x  axis  as  a.  Obviously, 

x  -  R  cos  a,  3  -  R  sin  a, 
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.v-v' 


where  i?  is  the  distance  from  the  center  of  the  area  to  the  point 
of  observation.  The  exponential  curve  in  (7)  will  now  be  written: 


exp  {iR  (*,  cos «  +  V k*  —  ftj  sin  a)}. 


Since  R  is  assumed  to  be  large,  when  k x  changes  the  exponential  curve 
will  be  an  extremely  rapidly  changing  function.  This  makes  it  possible 
to  use  either  the  stationary  phase  [2]  or  the  saddle-point  [3]  method 
to  calculate  the  integral.  In  this  case  the  principal  input  to  the 
integral  will  be  only  the  values  of  kx  lying  near  the  saddle  points 
kx  =  v,  which  are  determined  from  the  equation 

[(dldkj  (kj  cos  a  +  V &  —  k\  sin  =  0, 


from  whence 


v  "  k  cos  a  . 

To  calculate  the  integral  in  this  case,  in  the  first  approximation 
the  entire  subintegral  function,  except  the  exponential  curve,  may  be 
moved  beyond  the  integral  sign  when  kx  =  v;  the  expression  in  the  expo¬ 
nential  curve  may  be  expanded  in  a  series  according  to  the  powers  of 
kx  —v.  As  a  result,  we  obtain  from  (7): 


+®  . 

?=  2  S«  («)  <t>(k°  +  mp  —  k  cos  a)  e1**  x 

m«— co 

*  S  exp  i  —  (6 V - 1>)* M*'  (8) 

— CO  ' 


Here  Bm (kx,kz)  is  designated  by  Bm( a)  when  kx  =  k  cos  a  and  kz  «  k  sin  a; 
the  integral  in  this  case  is  reduced  to  the  well-known  Poisson  integral 
by  substituting  kx  —  v  =  e~ ^/4  •  and  it  yields 


Ak/R  sin  . 


In  view  of  the  fact  that  the  function  <t>(k0a  +  mp  —k  cos  a)  differs  from 
zero  only  in  the  narrow  interval  where  ix  «  p,  one  may  only  consider, 
of  all  the  terms  of  the  sum  in  (8),  only  that  term  for  which  m  is  equal 
to  a  whole  number  close  to  the  value: 
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(l/p)  (A:  cos  a  -  k°x ) . 


(9) 


As  a  result,  we  obtain: 


9  =  Bm  («)  sin  « <lJ(fcS  +  trip  —  k  cos  a)  exp  {—  (ir.fi)  -f  ikR)  V ~blfo  ( 10) 


Thus  the  diffraction  field  in  this  case  will  be  a  cylindrically 
diverging  wave  of  a  directional  nature  in  the  form  of  separate  lobes, 
the  maxima  of  which  lie  in  the  directions  of  a  =  a^,  which  are  determined 
by  the  equality: 


Asln<xm=  k*  +  mp,  m~ 0,  ±1,'  ±2,... 


(11) 


In  these  directions  the  argument  of  the  function  4>  become*  zero, 
and  the  function  itself  is  maximum.  The  width  of  the  lobes  is 
determined  by  the  magnitude  of  the  interval  Ax.  =  2r/D,  in  which 
the  function  is  notably  different  from  zero  [see  (4)  and  (5)]. 

In  conformity  with  this,  the  boundary  of  lobe' m  may  be  defined  as 
angle  a'm,  for  which 


kx  -}~  mp  —  k  cos  «„  sw  2n ID. 


(12) 
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From  (11)  and  (12)  we  obtain 


fcfcos  —  cos  a'm )  *  2-n/D; 

from  whence,  considering  the  closeness  of  the  values  of  angles  ct^  and 


m 


v.!m  =  \/D  sin 


(13) 


which  will  be  the  angular  half-width  of  each  lobe'.  The  number  of 
lobes  is  determined  by  the  interval  of  whole  values  of  m  where 
| cos  Opji  <  1  and  approximately  equal  to . 2A/1 .  The  results  obtained 
do  not  differ  in  principle  from  the  results  obtained  during  the  inves¬ 
tigation  of  diffraction  from  an  optical  grid  of  finite  dimensions  [4] . 

If  we  were  to  have  examined  diffraction  from  an  area  bounded  in  both 
directions,  we  should  have  obtained  results  analogous  to  those  presented 
in  [4]  (§  52)  for  a  two-dimensional  grid. 
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The  form  of  each-  lobe  is  determined  by  the  form  of  function  *£ 
in  (10),  and,  consequently,  according  to  (4),  the  form  of  function 
F(X.' .  In  the  case  of  a  limited  area  with  illumination  which  is  constant 
over  the  entire  area,  we  have 


« 

F(X)  - 1,  -£/2<X<£>/2, 
F(X)  =  0,  \X\>DI2.  . 


(14) 


and  from  (4)  we  find  that 


nrr-  OS) 


For  the  maximum  value  of  this  function  when' x  =  0  we  have  D/ 2n. 

The  function  also  has  auxiliary  maxima  which  may  be  disregarded  in  most 
cases.  The  value  of  the  function  in  the  first  auxiliary  maximum  when 
xD  =  3it  is  approximately  5  times  less  than  in  the  principal  maximum. 

The  half-width  of  the  principal  maximum  (from  the  center  of  the  maximum  to 
zero)  is,  according  to  (IS)', "equal  to 


;  Ax.  *  2 v/D  . 
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It  is  interesting  to  note  that  in  the  case  of  an  infinite  surface, 
as  is  evident  from  (14,  1)  and  (21,  1)  and  the  expressions  (32,  1)  and 
(33,  1)  analogous  thereto,  that,  in  the  case  of  electromagnetic  [waves], 
when  k™  0,  the  amplitude  of  the  corresponding  spectrum  increases 
infinitely.  In  the  case  of  a  finite  area,  the  amplitude  of  the  spectrum, 
in  accordance  with  (10),  is  proportional  to  B  (a)  sin  a,  and,  by  virtue 
of  the  fact  that  a  •+  0  when  k3  -*■  0,  the  amplitude  is  finite. 

The  above-obtained  results  are  valid  for  sufficiently  remote  points 
lying  within  the  Fraunhofer  zone  whore  the  condition 


R  »  kD2 


(16) 


is  satisfied.  This  condition  is  well  known  in  optics.  This  comes 
about  because  of  the  requirement  for  slowness  in  changing  the  function 
4>(fc0,.  <■  rirp  —  cos  a)  as  compared  to  the  exponential  curve  in  (8) ,  which 
is  necessary  in  order  for  it  to  be  possible"  to  move  this  function  beyond 
the  integral  sign. 
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2.  Diffraction  on  a  Sinusoidal  Surface 


When  examining  specific  examples  we  shall  limit  ourselves  to 
the  simplest  case  of  a  one-dimensional  irregularity,  assuming  that 
function  Z  depends  only  on  X.  In  the  case  of  an  infinite  surface, 
the  diffracted  wave  will  consist  of  the  aggregate  of  the  flat  waves, 
the  directions  of  distribution  of  which  will  be  set  by  the  components  of 
the  wave  vector  In  the  case  of  a  finite  area  of  width  D,  the 

diffracted  wave  wave  will  be  a  cylindrical  wave  with  a  characteristic 
of  directionality  consisting  of  the  aggregate  of  the  lobes,  the  maxima 
of  which  lie  in  directions  The  angular  half-wi'dth  of  the  lobes 

is  given  by  formula  (13).  The  amplitude  of  the  wave  in  the  center  of 
the  lobe,  according  to  (10)  and  (IS),  is  given  by: 


Am  —  [Djir:)  y-k/RBm(«m)  Sinam.  (17) 


In  the  acoustic  case  Bnrfa^)  is  found  from  (21,  1).  In  the 
electromagnetic  case,  when  the  horizontal  polarization  is  (Ey  ^  0, 

Ex  =  Ez  =0),  Bfh  (c^jj  is  calculated  frcm  (33,  1)  if  it  is  assumed  there 
that  q  *  0  at  the  same  time  that  Bjjj  =  =  0.  We  shall  note  that 

here  B %  (c^)  coincides  with  Bm  ( a mJ  for  the  acoustic  case  if  V  is  taken 

to  be  -1  in  (21,  1)  --  as  should  have  been  expected.  In  the  same  manner, 
in  the  case  of  vertical  polarization  {Hu  ^  0,  Hx  =  Hz  -  0) ,  for  a  magnetic 
field  we  shall  obtain  the  formula  (21,  1),  where  it  is  necessary  to  take 
V  to  be  equal  to  1.  However,  since  in  the  future  we  will  be  interested 

only  in  the  amplitudes  of  the  waves  and  not  in  their  phases,  the  cases 

of  V  =  1  and  V  -  -1  are  fully  equivalent. 

Thus  the  function 


f  "  \Bm  (S?)  I  sin  am  (18) 

characterizes  the  angular  dependency  of  the  amplitude  of  the  scattered 
waves  for  the  plane  problem  of  both  acoustic  and  electromagnetic  waves. 
Making  the  exchange  in  (21,  1) 


kg  1:3  k  Sin  C tmi  kjt  — ■  k  COS  0 tan 
kg'  —  —  Atsina„,  =*  k  cos «o 


(where  cto  is  the  slide  angle  of  the  incident  wave),  and  having  assumed 
that  |y)  =  1,  we  obtain: 
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(19) 


» .  •.  ^  cos  (am  *•)  i  DO  /„  \  t 

/(“»*)  =  "Tin  «’  +  sin  a7~ 1  Bm  m)  1  ’ . 


where,  according  to  (17,  1),  5°  (c^)  are  the  coefficients  of  expansion 
of  the  function:  m 


+00 


exp  {—  /^(slnaM  +  sin  a0)  Z (A")}  =  (*«)  e,"'pX’  (20) 


We  shall  switch  to  an  examination  of  scattering  from  a  sinusoidal 
surface:  in  this  case  we  have 


Z  =  a  cos  pX  . 


(21) 


As  is  known,  there  occurs  the  relationship 


+C0 


e—  *o»^x  »  ^  '(-<)"  (s)  e,mpX. 


(22) 


where  is  a  Bessel  function  of  order  m.  Therefore,  having  designated 


s  =  ka  (sin  am  ■(-  sin  «0) 


(23) 


and  compared  (20)  and  (22),  upon  calculating  (19)  we  obtain 


(am)  =s  ( (s). 


1 24) 


Substituting  (24)  into  (20),  we  obtain  the  final  formula  for 
wave  amplitude  in  the  maximum  leaves,  with  an  accu’  acy  of  up  to  the 
constant  factor: 


J  K  ml  sin  a  +  ssn  on  Jn  W- 


(25) 


Figures  1-3  present  graphs  of  the  directionality  charactf ristics  of 
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a  diffracted  wave  for  two  incident  wave  slide  angles  (45°  and  10°) 
when  there  are  various  values  for  6  =  2tt a/X,  where  a  is  the  amplitude 
of  the  sinuosity,  hhd  is  the  length  of  the  radiation  wave.  Tne  jergt.. 
of  the  wave  of  sinuosity  A  is  the  s arc  in  ail  cases  and  is  equal  tc 
A  =  10)..*  However,  the  width  of  the  lobes,  which  only  the  area  dimensi.r 
influence,  is  done  very  approximately  on  the  graphs.  A  view  of  the 
uneven  surface  is  depicted  in  the  lower  part  of  each  Figure.  Where 
the  amplitudes  of  sinuosity  are  very  low,  the  wavy  line  is  depicted 
as  a  straight  line.  The  arrows  on  the  graphs  indicate  the  directions 
of  the  incident  and  reflected  wave  where  specular  reflection  exists. 

The  force  of  the  sound  in  decibels  is  indicated  along  the  radius  in 
all  the  graphs. 

*The  length,  of  the  area  is  taken  as  equal  to  Uh: 


Figure  1:  6  *  2Tta/X  =  10;  A/^-  =  6.28 


We  see  that  even  at  relatively  low  sinuosity  (6  =  10,  i.e.,  a  = 

1.6.X,  Figure  1)  diffusion  scattering  almost  occurs.  Even  at  low 
amplitudes,  when  6  =  6  ar.d  6  =  3,  ir.nx.‘r.wun  radiation  scattering  occurs 
not  in  tne  direction  of  specular  reflection,  but  at  angles  equal  to 
100°  and  72°,  correspondingly.  This  is  easy  to  understand,  using  the 
idea  of  ra/s.  Rays  specularly  reflected  from  the  point  of  inflection  of  the 
sinusoid  go  in  these  directions.  In  view  of  the  fact  that  the  curvature 
of  the  surface  is  equal  to  zero  at  the  points  of  inflection,  the  rays, 
upon  reflection,  experience  the  least  separation,  as  a  consequence  of 
which  the  amplitude  of  the  corresponding  waves  is  maximum. 

The  situation  where  a  wavy  surface,  even  one  differing  very  little 
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from  flatness,  yields  a  noticeable  deviation  from  specular  reflection, 
as  is  evident,  for  example,  from  Figures  5  and  6,  is  interesting. 

Only  in  Figure  7,  where  the  amplitude  of  the  sinuosity  is  209  times 
less  than  the  length  of  the  radiation  wave  and  2090  times  less  than  A, 
may  the  reflection  be  considered  specular,  since  the  side  lobes  are 
approximately  30  db  less  than  the  principal  lobe. 

Where  the  slide  angle  of  the  incident  wave  is  10°,  the  scattering 
is  more  directional  (compare,  for  example,  Figures  3  and  9),  given  the 
same  6  values,  than  is  the  case  at  4Sg.  This  is  a  confirmation  of  the 
well-known  fact  that  at.  low  slide  angles  the  reflection  from  a::rough 
surface  is  more  directional. 

We  shall  also  note  that  in  all  cases  to  which  the  above-presented 
figures  apply,  the  theory  applicability  conditions  indicated  in  §  5,  1 
are  satisfied.  Specifically,  the  directions  of  the  lobes  depicted  in 
the  Figures  nowhere  greatly  differs  from  the  directions  of  possible 
geometric  reflections  of  the  rays. 


3.  A  Comparison  of  Sinusoidal  Surface  Solutions  Obtained  by  Various 
Methods 

In  the  case  of  a  sinusoidal  surface  one  may  obtain  a  solution  to 
the  problem  by  another  methods  Although  the  analysis  of  this  solution 
is  most  complicated,  it  is  useful  to  compare  it  to  our  solution  for 
the  case  where  the  amplitude  of  the  sinusoid  is  small  in  comparison 
to  wave  length. 

For  the  sake  of  simplicity  we  shall  examine  the  case  of  a  wave 
striking  the  surface  perpendicularly.  We  shall  again  write  the 
equation  for  the  surface  in  the  form 


Z  =  a  cos  pX  . 


We  shall  assume  the  surface  to  be  a  reflecting  surface,  assuming  that 
the  boundary  condition  =  0  is  fulfilled  thereon.  In  electrodynamics 
this  case  corresponds  absolutely  to  the  conducting  surface  and  to  the 
vector  E  which  is  directed  parallel  to  the  y  axis. 


We  shall  write  the  total  field  of  the  incident  and  diffracted  waves 
in  the  form2 


where 


?  =s  4.  AoC1**  +  cos  px  +  AtC'»*x  cos  2px  -f . . . , 


P„==J/£J  —  ntp1,  n  t=»  1,2,... 


(26) 


2  Qur  method  of  obtaining  a  precise  solution  is  essentially  a  generali¬ 
zation  of  Re^ley's  discussions  [5]. 
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This  kind  of  full  field  representation  coincides  fully  with 
the  separation  of  the  diffracted  field  into  spectra  [see,  for  example, 
formula  (14,  1)]  obtained  in  [1],  Here  each  term  of  the  formula  (26), 
beginning  with  the  third,  corresponds  to  two  spectra,  which  is  easy  to 
envision  if  the  cosines  are  expressed  by  exponential  functions. 

Substituting  the  surface  equation  Z  -  a  cos  pX  into  equation  (26) 
and  letting  9  =  0,  we  obtain: 


A0  +  exp  {-  2 iak  cospA'}  +  Ax  exp  {-  /s,  cos  pX)  cos  pX  + 

+  At  exp  {—  i$t  cospX)  cos  2 pX  + . . .  =  0,  (27) 

where 

(28) 


In  equation  (27)  the  exponential  • urves  may  be  replaced  by 
sums,  if  one  considers  that,  analogously  to  (22),  we  have 


c-i,  co,m  ^  ^  ^  ^  cos/J/y  ~  2 J3  (s)  cos  2 pX  +  . . . 

*  •  • 

In  the  expression  obtained  after  this  [operation]  the  products  of 
the  cosines  must  be  expressed  by  the  cosines  of  the  multiple  angles. 
It  is  also  necessary  to  group  the  terms  containing  the  cosines  of 
identical  arguments  as  factors.  Equating  the  co'efficients  before 
these  cosines  to  zero,  we  obtain  the  following  infinite  system  of 
equations  for  determining  the  coefficients  AQl  AXi  ...3 


where 


and  finally 

"T  m/i 


<0 

2  "fiiM A»  =  Cmi  tn  —  0,  1,  2.. . , 

rt-0 

(29) 

£0  1=5  Jo  (2tfA),  Too  =  1» 

(30) 

Cm  —  2  (  i)  Jm  (2ak),  nx  =  1 ,  2, ... 


(-  <r+"  [J*+m  (Sn)  +  (-  («„)],  m,  n  =  1,  2, . . ;  (31) 


Investigation  of  the  system  of  equations  (29)  is  an  extremely 
complicated  problem.  It  would  be  of  interest,  first  and  foremost, 
to  investigate  the  convergence  of  the  series  (26)  where  the  coefficients 
Aqi  A2,  ...  are  to  be  found  from  (29).  We  shall  limit  ourselves 

here  to  an  examination  of  the  case  of  small  irregularities  (ak  <■  <  1) , 

3  We  omit  here  computations  reduced  principally  to  a  transformation* of 
the  order  of  summing  in  compound  sums. 
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when  the  coefficients  decrease  very  rapidly  as  n  increases  [as  (ak)n], 
and  we  shall  compare  the  results  obtained  with  the  results  of  the 
above-presented  theory.  Retaining  in  the  expressions  for  An  only  the 
lowest  powers  of  ak,  we  find  from  (29)  that 


A$  =  —  1,  — 2i 'ak,  Ai  —  (a&)*£l  2a*}' 

A  - -t  W  [I  - 


(32) 


Changing  (26)  to  the  form  of  (14,  1)  by  means  of  substituting 
exponential  curves  for  cosines,  one  may  convince  oneself  that  Aq,  A\, 
Ai,  43,  ...  must  coincide  with  the  values  for  Sq,  2 Bj,  2B2,  2 B3,  ... 
The  latter  may  be  found  according  to  formulas  (26,  1)  and  (24).  Here 
it  is  necessary  to  take  into  account  the  fact  that  V  =  -1  with  our 
limited  conditions.  Besides  that,  since  the  incidence  is  assumed  to 
be  perpendicular,  =  -k,  kfl  -  0.  As  a  result,  we  obtain 


Bm~  ~{k  I  k?)  (—  i)m  Jm  [ak  (sin  am  +  sin  «#)] 


(33) 


and,  finally,  considering  the  smallness  of  ak  as  well  as  the  relation¬ 
ships  =  k  sin  am  and  aQ  =  rr/2 ,  we  find  that 


»  i 


Here 


2flj  =  aTiTHT,  ('T")?  o  +  sin  «>)*-. 


sin  «m  --£-  =  ]/ 1  ~ 


,  >n  =  1,  2,  3. 


Comparing  the  sequences  of  (32)  and  (34)  we  see  that  they  will 
coincide,  if  we  disregard  (mp)2  in  comparison  with  k2.  As  a  matter  of 
fact,  here,  on  one  hand,  sin  =  1,  and,  on  the  other  hand,  in  con¬ 
formity  with  (28),  8m  ;  0  and  we  obtain: 


A0  =  Ba  =  - 


I,  A,  =  2#i  =  2/n/v’, 


At  =  2  B3  =  (ak)-,  As  —  —  (i f/3 )  (ak)3 ...  ( 35) 

Thus  the  validity  of  our  solution  in  the  case  being  examined 
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is  limited  by  the  condition  that  (mp)2  «k2 ,  which  is  equivalent  to 
the  condition  that  sin  =  1,  i.e.,  =  tt/2 .  In  other  words,  we 

shall  obtain  correct  amplitudes  only  for  waves  scattered  in  directions 
close  to  perpendicular.  This  result  confirms  the  conviction  expressed 
in  §  5,1  above  thot  the  method  of  calculation  is  valid  only  for 
waves  scattered  in  directions  close  to  the  directions  of  geometri- 
reflection  from  various  parts  of  the  surface.  In  the  case  being  examined, 
in  view  of  the  smallness  of  the  amplitudes  of  the  sinuosity,  the  directions 
of  geometric  reflection  will  be  grouped  around  the  perpendicular  to  the 
average  surface  level. 


4.  Wave  Superposition 

We  shall  examine  the  case  where  a  wavy  surface  is  formed  by 
superpositioning  several  waves. 

We  shall  assume  that  we  have  solved  the  problem  of  diffraction 
on  a  surface,  the  form  of  which  is  given  by  a  certain  function  Z(X). 
This  means  that  the  coefficients  of  S°m  (o^)  have  been  found  in  an 
expansion  of  (20),  The  question  is  asked,  how  do  these  coefficients, 
and,  consequently,  the  amplitudes  of  the  diffracted  spectra  as  well, 
change  if  a  sinuosity  defined  by  the  cosine  curve  a  cos  ypX  and  having 
a  period  y  whole  number  of  times  less  than  the  period  of  function  Z(X) 
is  superimposed  on  this  surface? 

We  shall  designate  the  new  coefficients  which  are  analogous  to 
as  c°m  •  T^e  latter  ir,ay  be  found  from  an  expression 
which  is  analogous  to  (20) : 

exp  {—  ik  (sir.  xm  -f  sin  a0)  [Z  (AO  +  n  cos  vp.YJ)  =  ^ Cm  (^v)  e'mpx, 

m 

which,  when  (20)  and  (23)  are  computed,  is  written  in  the  form 
c~ht0"pX  +|  Blc,npX  =  g  C°me,npX. 

(TIU»— CO 

Substituting  the  series  (22)  for  the  exponential  curves  in 
the  left-hand  portion,  multiplying  out  the  series  and  grouping 
terras  correspondingly,  we  obtain 


2  **  2  <- 
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*«*£*'«<•*  ’Jifcvw  r-^c-o  >vv»>vv^*«.>  v- 


from  whence 


C2.M-  S  (-0V*(s)  *•-*,(«„). 


Figure  14  depicts  the  scatter  field  in  the  case  where  the 
sinuosity  is  a  superpositioning  of  two  sinusoids.  Calculation  is 
accomplished  with  the  help  of  formula  (36).  The  angle  of  wave 
incidence  is  taken  to  be  45°.  Besides  that,  it  is  taken  that 


Sj  =  ka^  =  3,  and  62  =  ka2  =  2, 

where  and  a2  are  the  amplitudes  of  the  sinusoid.  As  in  the 
preceding  cases,  it  is  assumed  that  A/X  =  10. 


angles  in 
degrees 
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Figure  lk:  <5  =  3,  <5j  “  2 


Figure  15:  individual  sinusoids  comprising  the  surface  shown 
in  Figure  \k 
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A  view  of  the  scattering  surface  is  depicted  in  the  lower 
portion  of  the  Figure  [14].  Individual  sinusoids,  the  superposition 
of  which  results  in  this  surface,  are  depicted  in  Figure  15. 
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5.  Diffraction  on  a  surface  which,  [when  viewed]  cross-sectionaiiy , 
is  a  broken  line 

We  shall  examine  diffraction  on  a  surface  defined  by  the 
equations  (Figure  16) : 


ft 

& 


2  =  6(1  -H*/A),  -  A  /  2  <  A!  <  0, 

2  =  a(l~4X/A),  0 ^ X A / 2. 


(37) 


U&  * 

I  S' 

I  | 


For  the  expansion  coefficients  of  (20)  we  have,  according  to  the 
known  formulas  for  Fourier  series  coefficients, 

,  A/2 

Bl  («„,)  =  ™  S  exP  ik  (sin  **  +  sin  «o)  Z{X)~  irnpX]  dX  (38) 


Figure  (6:  Cross-sectional  view  of  a  "broken- I i ne"  surface 
or,  taking  into  account  (37)  and  (23), 


b;,  («n)  s=s  ^  exp  |  —  is(l  +  —)  ~  impx}  dX  + 

-A/2 

'  A/3  . 

+  ™  $  exP  {  “  * (J  “x)“ top*}**- 

$ 


Having  compJU’ted  the  integration,  we  obtain 

*  («M)  =  sin  (x  +  s). 

V mj ) ~s 


(39) 
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where  the  relationship  p  =2n/A  is  also  taken  into  account. 


Figure  17:  Scattered  wave  directionality  characteristic 
according  to  formula  (40) 

6  =  2na/A,  A/a  =  20.9 


Having  substituted  (39)  into  (19),  we  find  for  the  amplitudes 
of  the  waves  at  the  lobe  maxima  with  an  accuracy  up  to  a  constant 
factor  that: 


• //„  «*fI-cos(aw-a#)l  .  fmn  .  \ 

/{m) - 7*  sln  (t  ' +  5) ■ 


(40) 


Figure  17  depicts  the  scattered  wave  directionality  characteristic 
calculated  according  to  formula  (40) .  The  angle  of  incidence  is  taken 
to  be  45°,  A/A  =  10,  and  <5  -  2rra/A  =  3. 

As  may  le  seen,  the  scattered  wave  has  a  character  close  to  that 
of  diffusion.  The  lobe  directed  at  an  angle  of  approximately  66°  has 
maximum  amplitude.  This  lobe  corresponds  to  specular  reflection  from 
the  flat  areas  forming  the  left  slopes  of  each  peak  in  the  uneven 
surface  (for  example,  the  sector  —  A/2  X  ^  0  in  Figure  16). 

The  lobe  directed  at  an  angle  of  approximately  23°  and  corresponding 
to  reflection  from  the  right  slopes  is  also  most  intensive. 

As  was  already  indicated  above  (see  §  5,  1),  our  theory  is  able 
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to  give  correct  values  for  lobe  magnitudes  only  for  angles  which 
do  not  differ  too  greatly  from  these  two  directions  of  geometric 
reflection.  At  angles  greater  than  90°  j.t  may  not  be  expected  that 
the  dimensions  of  the  lobes  will  be  correct.  They  are  depicted  in 
Figure  17  only  to  illustrate  the  resulis  of  the  theory,  which  is 
based  on  Kirchhoff’s  principle. 

The  scatter  diagrams  presented  here  were  constructed  by  Ye. 
Zezyukina;  she  also  did  all  the  necessary  calculations.  In  connection 
with  this,  I  express  my  gratitude  to  her. 
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